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Abstract: The textbook treatment in that the wave function of a dynamical 
system is expanded in an eigenfunction series is investigated. With help of 
an elementary example and some mathematical theorems, it is revealed that 
in terms of solving the time-dependent Schrodinger equation the treatment 
involves in divergence trouble . The root reason behind the trouble is finally 
analyzed. 
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Solutions of the time-dependent Schrodinger equation, formal ones and 
numerical ones, are of considerable interest in the past and in the present. 
It is recognized that they are related to many interesting and essential dy- 
namical processes in the quantum realm. 

Being aware of difficulties with solving the time-dependent Schrodinger 
equation directly, the textbook treatment takes the detour-type approach 
in that a wave function is expanded in an eigenfunction series and thus the 
dynamics of coefficients of the series, instead of that of the wave function, 
becomes the focus of the theory Q. The procedure was put forward ini- 
tially by Dirac at the early time of quantum mechanics and has since been 
accepted as one of standard methods in many perturbative and nonpertur- 
bative approaches. 

Though not enough attention was received, discussions and suggestions 
concerning the Dirac approach existed for decades H||. I n some of our 
recent papers, it was argued that even for today's physicists there were 
still things worth careful consideration j|]. The purpose of this paper is to 
examine the pure mathematical aspect of the subject. By the term "the pure 
mathematical aspect" we mean that (i) the foundation of solving the time- 
dependent Schrodinger equation, not any concrete perturbative procedure, 
will be of interest; (ii) only simple examples that can be investigated both 
numerically and analytically will be discussed; (iii) except the energy concept 
no physical arguments, particularly those of gauge, will be related. 

The basic assumption of the standard treatment is that a system's wave- 
function can be expanded in a series of the following form 

VP(i,r) =^C n (t)e-^/^ n (r), (1) 

where C n (t) stands for a coefficient that is a constant before the initial 
time and then becomes time-dependent after it and \£n( r ) is one of the 
eigenfunctions satisfying the eigenvalue equation 

H V n {r) = E n V n (r), (2) 

in which Hq, according to the convention, represents the Hamiltonian of the 
system at the initial time. Provided that all eigenfunctions are adequately 
normalized, the coefficients of the expansion must obey the normalization 
condition 

oo 
£|Cn(*)| 2 = l. (3) 



Inserting (|l|) into the time-dependent Schrodinger equation 

ih— = H(t,r)V, (4) 

multiplying both the sides by e "'' \&*(r) and integrating the resultant 
equation term by term, we arrive at a set of coupled ordinary-differential 
equations 

ifr—fT = E ClVnl ex P( iuJ mt), (5) 

where co n i = (E n — Ei)/fi and V n i stands for the matrix element of the 
Hamiltonian variation V{t) = H(t) — Hq. It is almost unanimously be- 
lieved that the equation set (|5|) is completely and exactly equivalent to the 
Schrodinger equation (Q) and, at the same time, enjoys some advantages. 
One of the advantages is that the equation set appears to be numerically 
solvable and some efforts are indeed in this direction ||. According to the 
customary thought, there is nothing to be discussed concerning the validity 
of (H) because (i) no approximation has been introduced; (ii) no perturba- 
tive procedure has been invoked; (hi) no physical interpretation has been 
imposed. 

Nevertheless, we will, in this paper, show that the formalism outlined 
above is not truly reliable since the basic expansion ([[]) is not a mathemat- 
ically well-behaved series. 

Before discussing the issue in an abstract way, let's apply the standard 
formalism to a one-dimensional oscillator and see where it leads us to. The 
oscillator is assumed to have the Hamiltonian 

where 

f 1 (t < 0) 

S(t) = I 1 + rft (0 < t < T) (7) 

[l + 7]T (t>T). 

(Several other types of time-dependent Hamiltonians were also tried. The 
situations were all alike as far as the subject of this paper was concerned.) 
It is well-known that the system expressed by (g) is a simple harmonic 
oscillator before t = 0. The quantum state before t = is expressed by 



E 

n=0 



C n0 e-^ n+1 ^ t N n (a)H n (ax)e-2 Q x , (8) 



where u = ^/kjm, a = (mk/h 2 ) 1 / 4 , N n = [a/(^2 n n\)] 1 / 2 and H n is the 
nth Hermite polynomial, which is governed by the iteration relation 

Hn+iiO = 2£H n (0 - 2n J ff„_ 1 (0 (9) 

with 

H (0 = 1, ffi(0=2C, H 2 (Z) = 4f-2 (10) 

We will assume that the system is in the ground state before t = 0, 
namely we have 

Co(0) = l and C n (0) = (n ^ 0), (11) 

and try to calculate the quantum state at the later times. 
To make the numerical work as simple as possible, we set 

77 = 1, T = l, m = l, k = l, h = l. (12) 

Under these conditions, the Hamiltonian variation V(t) becomes 

x 2 t/2 (0 < t < T) 
x 2 /2 {t > T). 

By virtue of the iteration relation (|9|) , we obtain the following relations 

(2j|x 2 /2|2j) = (j + 0.25) (14) 

and 

(2j\x 2 /2\2j + 2) = 0.5 x [(j + 0.5) (j + 1)]3. (15) 

With the notation C n = C r n + iC l n , the equation set (||) becomes 

dC r 

-£L = [C l 2j _ 2 cos(2t) + C r 2j _ 2 sm(2t)]V 2j „ 2 , 2j 

+C l 2j V 2j>2j + [C l 2j+2 cos(2t) - C r 2j+2 sm(2t))V 2j;2j+2 



V(t) = { 2 L ).. ^ ' (13) 



and 



dCii 



iU ~{[C r 2] - 2 cos(2t) - C\ 3 _ 2 sin(2i)]F 2i _ 2 , 2j (1?) 

+C r 2] V 2j , 2j + [C r 2j+2 cos(2t) + Ci j+2 sin(2t)]V 2j , 2j+2 }, 



where j runs over all non-negative integers. Note that in the equations 
above any quantity having a negative subindex actually vanishes. Namely, 
C r _ 2 = Ci 2 = and y_ 2 ,o = 0. 

In numerical mathematics the Runge-Kutta algorithm || says that for a 
vector-type differential equation 

y = f(y,i), (is) 

we can construct the solution according to the formula 

yi+i = yj + g( z i + 2z 2 + 2z 3 + z 4 ), (19) 

where 

yo = y(*o), yi = y(h), •••,yj = y(tj),--, (20) 

with tj — tj— 1 = h, and 

zi = hf(yj,tj) 

z 2 = M(y j +z 1 /2,t j + h/2) 
z 3 = M(y j +z 2 /2,t j + h/2) 
z 4 = M(yj +z 3 ,tj + h). 

Although the algorithm proves very effective in classical mechanics ||, it 
encounters difficulties one by one in dealing with ( |l6| ) and ([l?]). Firstly, it is 
found that the number of involved variables, namely nonzero C 2 a and CL, 
increases rapidly. At t = 0, only Cq = 1 is there. At t = nh, we have (8n + 
2) nonzero variables. Secondly, when the code, adapted to the increasing 
number of variables, is run in a computer, it is found that if the Hamiltonian 
V(t) is relatively large, the normalization condition J2 |Cn(^)| 2 = 1 will soon 
be broken down. Figure 1 illustrates that the breakdown takes place within 
an oscillation cycle in respect to the perturbation described by ( |i~3| ) and 
h = 0.001. Thirdly, it is noted that whenever the normalization condition 
runs into difficulty, the average energy of the system 

(E)=Y / E n \C n (t)\ 2 (21) 

tends to infinity in a similar way as shown in Figure 2. (The code is available 
if more details are interested.) 

If we study the numerical calculation in an analytical manner, more 
puzzling things surface impressively. Equations (|0]) and (|l7|) suggest that 
if the system is in the ground state initially, it will, at the next step, be 



partly in the second excited state owning to the coupling of the coefficients. 
At the next next step, it will be partly in the fourth excited state. Since 
the Hamiltonian variation V(t) (O) has a finite value as t — > oo, the process 
will never terminate and the state of the system will vary forever. However, 
elementary quantum mechanics |/J asserts that the system will be settled in 
a stationary state after t = T. The contrast becomes more striking if we 
look at Figure 2, in which the "energy" curve rises and falls at the times 
t > T = 1 disregarding any theory about stationary states. 

Another perplexing thing is related to energy more directly. According 



to the formulas (||) and (21), a ground-state system cannot lower its energy 
under any circumstances. The physical intuition, however, tells us a different 
story. For instance, if the Hamiltonian variation of our oscillator takes the 
following form 

V(t) = -xH/2 (t > 0), (22) 

the system, in contrast to the system affected by (|l3|) , should lose its energy 
since, roughly speaking, the perturbation ( p2[ ) makes the potential energy 
of the system smaller while leaving the kinetic energy unchanged. 

We now try to work out the puzzle related to the title of this paper: 
Why does the standard formalism, derived from the Schrodinger equation 
in a seemingly rigorous way, suffer from so many difficult things? For our 
purposes, the analysis will be done almost entirely from mathematical con- 
sideration. 

It should be pointed out that 

e -tfotA$ ( r ) ) e - iElt / h ^i{r), • • • , e-"^*/ R tf n (r), • • • (23) 

do not constitute an orthogonal-and-complete basis for a function /(£, r). 
Consider our one-dimensional case again. Many functions, such as 

t, (t + x), cos{tx)e-^ a2x \ e tx e~-2 a2x \ , (24) 

cannot be expanded in a series of the form (JSj) . That is to say, the practice of 
expanding the solution of a Schrodinger equation in a series of eigenfunctions 
is originated from an intuitive physical idea rather than from a stringent 
mathematical analysis. 

The major argument, which we wish to propose in this paper, is the 
following. Even if the wave function at any fixed time can be expanded in 
an eigenfunction series of the form 

X;C tt (t)tf n (r), (25) 



substitution of this expansion, or the equivalent series (JTJ) , into the time- 
dependent Schrodinger equation is still a questionable procedure. As a mat- 
ter of fact, there indeed exist mathematical theorems restricting the use of 
infinite series []|]. One of them states that only a series converging uniformly 
can safely be integrated term by term. Another of them indicates that to 
differentiate a series term by term the resultant series has to be uniformly 
convergent. An inspection can tell us that any known eigenfunction series, 
if given as a truly infinite one, behaves poorly in this respect and cannot 
be safely used in the time-dependent Schrodinger equation. To see this in 
an example, apply the standard treatment to the series (|8|). Assume that in 
the series 

C„o = —. — — r, such that V \C n0 \ 2 = —~ V - — -—^ = 1. (26) 

By taking the time derivative, the left side of the Schrodinger equation 
becomes 



ih—— = e 2 

dt it 



oo 

'e-^ n+l ^ wt N n (a)H n (ax) 



^n+1/2 



n + 1 






(27) 



For this step and the next step to make sense, in which the expression has 
been differentiated and will be integrated term by term, it is necessary to 
ensure the uniform convergence of the series included in the square bracket 
above. (Note that the common factor e~^ ax ' > 2 is not relevant to the treat- 
ment term by term.) By setting t = 0, this series becomes 



V ^±±llN n (a)H n {ax) ~ V /^f^" ( 28 ) 

^ n + 1 v ' y ' ^ (2 n n\) 1 / 2 

whose value rises very fast for large x (faster than e 2ax ). A simple analysis 
tells us that the series cannot be uniformly convergent since it has no upper 
limit on its domain. 

At this stage, it is desired to have explanation for the fact that the series 
form 

C n0 e- lE ^H n (v), (29) 

where C n o are pure time-independent constants, has safely been used for 
almost all stationary states. We believe that there are two facts responsible 
for it. The first one is that for practical cases the initial expansion of ( |29| ) 
takes the form 

N 

]TC n0 *n(r), (30) 

n=0 



£• 



where N is a finite number. The finiteness of the series is due to the fact 
that high-energy eigenfunctions are related to the wave function in the very 
remote regions, which are usually of no interest for the problem, in which 
a fixed accuracy is set up at the very beginning. The second one is that 
each term in ( p0| ) evolves in the Schrodinger equation independently (no 
coupling), and therefore the series, if finite initially, keeps finite in all treat- 
ments. Keeping these facts in mind, one may say that a stationary wave 
function can indeed be expanded in a series of eigenfunctions (with limita- 
tions still, see Ref. 4). 

The discussion above has almost already manifested that things will be 
different for nonstationary cases. We may still express the dynamical wave 
function as a finite series at the initial time, but the series soon becomes a 
truly infinite one. The strong coupling, as (||) exhibits, tends to make each 
of the coefficients have a significant value. This is disastrous in view of that 
any truly infinite eigenfunction series suffers seriously in the time-dependent 
Schrodinger equation, as has been unveiled. 

From a slightly different viewpoint, we can also see that the expansion 
(|l|) cannot work smoothly. It is obvious that after the initial time the eigen- 
functions defined by Hq are no longer adequate eigenfunctions and the eigen- 
frequencies defined by u) n = E n /h are no longer adequate eigenfrequencies. 
If we forcefully express the dynamical wave function in a series constructed 
from these "out-of-date" eigenfunctions and "out-of-date" eigenfrequencies, 
the coefficients of the series have to adjust themselves violently. Since the 
large scale (the entire space) and the fast variation (represented by e ~ mujt ) 
get involved, the adjustment has to be so violent that the normalization 
condition ^2 \C n \ 2 = 1 will definitely break down. 

We proposed two alternatives to the standard treatment Hfl(|. One is 
based on a perturbative approach and the other is based on a nonperturba- 
tive approach. In the nonperturbative approach, with help of "intermediate" 
eigenfunctions and eigenfrequencies a general procedure of solving the time- 
dependent Schrodinger equation is constructed. 
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Ministry, P.R. China. 
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Figure captions 



Fig. 1, The time behavior of the coefficients in terms of the normalization 
condition. Note that the period of the oscillator is 6.28 since u = 1. 

Fig. 2, The time behavior of the coefficients in terms of the average energy. 
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Figure 2 
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